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Abstract 

Background: Detection of molecular tumor heterogeneity has become of paramount importance with the advent of 
targeted therapies. Analysis for detection should be comprehensive, timely and based on routinely available tumor samples. 

Aim: To evaluate the diagnostic potential of targeted multigene next-generation sequencing (TIVl-NGS) in characterizing 
gastrointestinal cancer molecular heterogeneity. 

Methods: 35 gastrointestinal tract tumors, five of each intestinal type gastric carcinomas, pancreatic ductal 
adenocarcinomas, pancreatic intraductal papillary mucinous neoplasms, ampulla of Vater carcinomas, hepatocellular 
carcinomas, cholangiocarcinomas, pancreatic solid pseudopapillary tumors were assessed for mutations in 46 cancer- 
associated genes, using Ion Torrent semiconductor-based TM-NGS. One ampulla of Vater carcinoma cell line and one 
hepatic carcinosarcoma served to assess assay sensitivity. TP53, PIK3CA, KRAS, and BRAF mutations were validated by 
conventional Sanger sequencing. 

Resu/ts: JM-NGS yielded overlapping results on matched fresh-frozen and formalin-fixed paraffin-embedded (FFPE) tissues, 
with a mutation detection limit of 1% for fresh-frozen high molecular weight DNA and 2% for FFPE partially degraded DNA. 
At least one somatic mutation was observed in all tumors tested; multiple alterations were detected in 20/35 (57%) tumors. 
Seven cancers displayed significant differences in allelic frequencies for distinct mutations, indicating the presence of 
intratumor molecular heterogeneity; this was confirmed on selected samples by immunohistochemistry of p53 and Smad4, 
showing concordance with mutational analysis. 

Conclusions: IM-NGS is able to detect and quantitate multiple gene alterations from limited amounts of DNA, moving one 
step closer to a next-generation histopathologic diagnosis that integrates morphologic, immunophenotypic, and multigene 
mutational analysis on routinely processed tissues, essential for personalized cancer therapy. 
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Introduction 

Cancer inter-tumor and intra-tumor heterogeneity, a well-lcnown 
fact described by pathologists in the classification of tumors over the 
last two centuries, has finally risen to the forefront of clinical interest. 
Cancer genomics and transcriptomics studies have shown that 
tumors belonging to the same histotype display remarkable 
differences in their genetic assets; such inter-tumor heterogeneity 



is the basis of molecular subclassification witli clinical impact for 
targeted therapeutic approaches. It has also become clear that 
phenotypicaUy and genetically diverse clones of neoplastic cells may 
be juxtaposed within the same tumor[l ,2] . These clones are thought 
to be players in a branching clonal evolution scenario leading to the 
formation of metastases that are more aggressive and resistant to 
treatments than the primary tumor [1]. 
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The histological and immunohistochemical characterization of 
multiple samples from the same tumor can highlight the presence 
of subpopulations of neoplastic cells displaying peculiar morpho- 
logical and immunophenotypical features; this morpho-phenotyp- 
ical analysis of intratumor heterogeneity finds its natural comple- 
ment in a comprehensive characterization of molecular lesions 
within a cancer specimen. The sum of these data offers essential 
information to diagnose and subclassify cancers for the scope of 
determining prognosis and selecting tailored treatments [3] . 

The sequencing analysis of hotspot mutations in cancer-related 
g(;nes has tiuis Ixxome a usc'ful tool in selecting personalized 
therapy for many malignancies [4,5]. However, the use of 
conventional techniques for a wide molecular characterization of 
tumors is hampered by the high costs and time needed to assess 
multiple molecular alterations, and by the limited amount of tissue 
consisting in formalin-fixed paraffin-embedded (FFPE) biopsies 
and/or fine needle aspiration cytology. This calls for the 
implementation of companion diagnostic methods for (i) simulta- 
neously testing multiple genetic alterations and (ii) quantifying the 
molecular subclones, i.e. the amount of cancer cells harboring any 
different mutation. 

Massive parallel sequencing, also known as next-generation 
sequencing (NGS), has recently been introduced and is the most 
sensitive approach to index multiple genes starting from a limited 
amount of DNA [6,7]. In the present study, we assayed a targeted 
multigene NGS (TM-NGS) test in 35 FFPE samples from diverse 
upper gastrointestinal tract tumors to define its diagnostic potential 
in characterizing cancer molecular heterogeneity. 

Materials and Methods 

Tumor samples 

A series of 35 formalin-fixed paraffin-embedded (FFPE) samples 
from surgically resected neoplasms, representative of diverse upper 
gastrointestinal and hepatobiliopancreatic cancer t\'pes (Table 
SI), were assayed for intragenic mutations in 46 cancer-related 
genes by TM-NGS. The series included 5 intestinal type gastric 
carcinomas (GC), 5 pancreatic ductal adenocarcinomas (PDAC), 5 
pancreatic intraductal papillary mucinous neoplasms (IPMN), 5 
ampulla of Vater carcinomas (AVC), 5 hepatocellular carcinomas 
(HCC), 5 intrahepatic cholangiocarcinomas (ICC), and 5 pancre- 
atic solid pseudopapUlary tumors (SPT). These latter had matched 
fresh-frozen and FFPE samples available and served to assess the 
performance of TM-NGS on DNA from both sources. In addition, 
DNA from cancer cell line AVCl [8] and one hepatic 
carcinosarcoma [9], served to assess sensitivity of the TM-NGS 
mutational assay. 

Ethics 

A total of 35 samples from 35 patients, acquired by the ARC- 
Net biobank at the University and Hospital Trust of Verona - 
Italy, were used in the present study. Of these, the materials from 9 
patients have been collected with the written informed consent for 
their use in research under Program 1885 (creation of a biobank) 
protocol 52438 approved by the local ethics committee of the 
Integrated Unversity Hospital Trust of Verona on November 23rd 
2010. This approval covers biological material collection for the 
ARC-Net coordinated biobank of samples from all cancer patients, 
including neoplastic and associated local and distant normal tissue. 
Protocols for collection included informed consent, approved 
under this program, from the patient to collect residual tissue 
samples for molecular research. The program includes approved 
amendments to address the later regulatory issues of sensitive data 
in genomic studies and a separate informed consent for access to 



sensitive data. These informed consents, received from patients, 
are registered in the biobank database together with samples 
collected. Samples from the remaining 26 patients had been 
collected prior to November 23rd 2010. These samples were 
acquired by the biobank according to a protocol approved by the 
local ethics committee for use in research of residual pathological 
tissue under an amendment to the above mentioned program 
1885 (creation of a biobank). The protocol indicates the procedure 
to acquire and register these samples, and anonymize patient 
information where it was not reasonable to reconsent patients 
directly. Thus, the consent for use of samples from these 26 
patients was waived by the ethics commitee. 

DNA extraction and qualification 

Neoplastic cellularity was assessed by microscopic examination 
and, when below 50%, enriched by manually microdissecting four 
consecutive 10 nm thick sections. AH samples were microdissected 
excluding the five SPT, HCCl to HCC4, and ICCl. 

Genomic DNA from frozen or FFPE tissues was extracted using 
the QiAamp DNA Mini Kit or QIAamp DNA FFPE Tissue Kit 
(Qiagen), respectively. Purified DNA was quantified and its quality 
assessed using Qubit (Invitrogen Life Technologies) and Nano- 
Drop (Invitrogen Life Technologies) platforms [10]. DNA 
suitability for PGR downstream apphcations was further evaluated 
tiirough BIOMED 2 PGR multiplex protocol and tiie PGR 
products were evaluated by DNA 1000 Assay (Invitrogen Life 
Technologies) on the Agilent 2100 Bioanalyzer on-chip electro- 
phoresis (Agilent Technologies) [11]. 

Targeted Multigene Next Generation Sequencing of 
iVlultiplex PGR Amplicons 

Twenty ng of DNA were used for multiplex PGR amplification 
using the Ion AmpliSeq Cancer Panel (Life Technologies) that 
explores selected regions of the following 46 cancer-associated 
genes: ABU, ALK, AKTl, APC, ATM, BRAF, CDHl, 
CDKN2A, CSFIR, CTNNBl, EGFR, ERBB2, ERBB4, 
FBXW7, FGFRl, FGFR2, FGFR3, FLT3, GNAS, HNFIA, 
HRAS, IDH1,JAK2,JAK3, KDR/VEGFR2, KIT, KRAS, MET, 
MLHl, MPL, NOTCHl, NPMl, NRAS, PDGFRA, PIK3CA, 
PTEN, PTPNll, RET, RBI, SMAD4, SMARCBl, SMO, SRC, 
STKll, TP53, VHL. 

The quality of the obtained library was evaluated by the Agilent 
2100 Bioanalyzer on-chip electrophoresis (Agilent Technologies). 
Emulsion PGR was performed either manually or with the 
OneTouch DL system (life Technologies). Sequencing was run on 
the Ion Torrent Personal Genome Machine (PGM, Life Technol- 
ogies) loaded with a 3 1 6 chip as per manufacturer's protocol. Data 
analysis, including alignment to the hgl9 human reference genome 
as well as variant calling and filtering, was done using the Torrent 
Suite Software v.3.6 (life Technologies) with default options for the 
Ion AmpliSeq Cancer Panel. Variants were annotated using the 
SnpEfF software v. 3. 4 [12] and the NCBI mRNA Reference 
Sequences listed in Table S2. Alignments were visually verified 
with the Integrative Genomics Viewer; IGV v.2.3 [13]. 

DNA Sanger Sequencing and Immunohistochemistry 

To validate the mutations detected by TM-NGS, TP53 (exons 
5, 6, 7, 8), PIK3CA (exon 10), KRAS (exon 2) and BRAF (exonl 1 
and exon 15) specific PGR fragments were analyzed by conven- 
tional Sanger sequencing, as described previously [14—16]. The 
immunohistochemical expression of p53 (clone DO-1, Immuno- 
tech, dilution 1:50) and Smad4 (clone B8, Santa Cruz Biotech- 
nology, dilution 1:200) was tested as a surrogate validation of the 
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TM-NGS results for these genes. The protocol included depar- 
afEnation of 4- (im FFPE sections in xylene, rehydration via 
decreasing concentrations of alcohol down to pure water, non- 
enzymatic antigen retrieval in citrate buffer (pH 6.0) for ,30 min- 
utes at 95°C. Immunolabeling was developed u.sing the Novolink 
polymer detection kit (Leica Microsy.stems) according to the 
manufacturer instructions. Appropriate positive and negative 
controls were run concurrendy. 

Results 

TM-NGS yields overlapping results on DNA from frozen 
and paraffin samples and quantitates the mutated alleles 

We used 5 SPT for which matched fresh-frozen and FFPE 
samples were available to test the proficiency of the assay using 
DNA from FFPE tissues. This tumor type is ideal for this purpose 
because it is characterized by a monotonous composition with a 
neoplastic cellularity of about 70-80% and has a molecular 
hallmark consisting of a heterozygous CTNNBl mutation [17]. 
Both fresh-frozen and FFPE specimens were assessed for 
neoplastic cellularity by two independent pathologists and the 
extracted DNA was subjected to TM-NGS. The quantity of 
sequences obtained was similar for fresh-frozen and FFPE derived 
DNA (Figure 1). The CTNNBl gene mutation was detected in 
all samples; moreover, the aUehc frequency of CTNNBl mutation 
was consistent with the percentage of tumor cells as scored by the 
pathologists (Table 1, figure 1). 

TM-NGS is highly sensitive on DNA from both frozen and 

paraffin tissue 

The limit of mutation detection of TM-NGS on DNA from 
fresh-frozen samples was assessed using DNA from AVC 1 cancer 
cell line with known mutations [8] serially diluted with non-tumor 



DNA from a commercial source (Universal unmethylated DNA, 
Chemicon Int., BUlerica, MA). AVCl cell line harbors the 
following homozygous variants: KRAS G12A and CTNNBl 
S45F somatic mutations and the nonpathogenic TP53 P72R 
variant. The commercial DNA is heterozygous for the common 
TP 53 P72R nonpathogenic polymorphism. AVCl and commer- 
cial DNA were mixed to obtain samples with a decreasing relative 
AVCl DNA content: 50%, 25%, 20%, 15%, 10%, 7.5%, 5%, 
2.5%, 1% and 0%. Twenty ng of each dilution point were 
subjected to TM-NGS with the Ion AmphSeq Cancer Panel. The 
three variants were identified in all samples containing AVCl 
DNA down to a frequency of 1 % (Figure 2A). The ratio between 
tumor DNA c:ontent and allelic frequency was approaching one 
for CTNNBl and TP53 variants, while KRAS mutation showed a 
higher ratio (2.03±0.18). This is consistent with the previous 
AVCl characterization showing copy number gain of chromo- 
some 12p, where KRAS resides [8]. 

To assess whether TM-NGS has the same detection limit on 
DNA from FFPE, we used two different tumor components from a 
previously characterized hepatic carcinosarcoma [9]. DNA from 
the microdissected hepatocarcinoma and sarcoma components 
was mixed to obtain samples with a decreasing relative hepato- 
carcinoma DNA content: 100%, 95%, 90%, 75%, 50%, 25%, 
10% and 0% (Figure 2B). Twenty ng of each dilution point were 
subjected to TM-NGS with the Ion AmphSeq Cancer Panel. 
Three known different genetic variants of the hepatocarcinoma 
component were used to assess the assay sensitivity: ABLl intronic 
g.l64164 G>T (variant frequency hepatocarcinoma = 24% sar- 
coma =0%), P1K3CA H1047R (variant frequency hepatocarci- 
noma = 28% sarcoma = 0%) and TP53 F109C (variant frequency 
hepatocarcinoma = 87%; sarcoma = 55%). The three mutations 
were identified in all samples containing hepatocarcinoma DNA 
down to a frequency of 2%, corresponding to the frequency of 
ABLl gene mutation in the 10% diluted sample. Moreover, the 
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Figure 1. Targeted multigene-next generation sequencing analysis of five solid pseudopapillary tumors. Depth of sequencing 
(coverage) of targeted regions analyzed in 5 matched fresh-frozen and formalin-fixed paraffin embedded samples of solid pseudopapillary tumor. 
Dots describe the coverage of each target sequence per sample; the quantity of sequences obtained was similar for fresh-frozen (F) and formalin- 
fixed paraffin embedded (P) derived DNA. 
doi:1 0.1 371/journal.pone.01 04979.g001 
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Table 1. Concordance between tumor cellularlty and CryVA/6//p-catenin mutation prevalence detected by deep sequencing In five 
solid pseudopapillary tumors (SPT). 





Case ID 


Sample type 


Tumor cells 


CTNNBl mutation allelic frequency 


SPTl 


Frozen 


80% 


39% 


SPT1 


FFPE* 


80% 


39% 


SPT2 


Frozen 


80% 


42% 


SPT2 


FFPE 


80% 


42% 


SPT3 


Frozen 


80% 


46% 


SPT3 


Frozen 


80% 


39% 


SPT4 


Frozen 


70% 


27% 


SPT4 


Frozen 


70% 


42% 


SPT5 


Frozen 


70% 


28% 


SPT5 


FFPE 


70% 


36% 



*FFPE = formalin-fixed, paraffin embedded specimen. 
doi:l 0.1 371/journal.pone.Ol 04979.t001 



detected variant frequency was consistent with the expected value 
(computed from the mutation frequency in each component and 
the percentage of hepatocarcinoma and sarcoma components at 
each point) at hnear regression analysis, showing that this assay 
may quantitate the actual allelic frequency of a somatic mutation 
in a given FFPE sample. 

TM-NGS describes intertumoral and intratumoral 
molecular heterogeneity 

We applied the Ion AmpliSeq Cancer Panel to a series of FFPE 
samples from 30 additional upper gastrointestinal tract tumors; the 
series consisted of 5 GC, 5 PDAC, 5 IPMN, 5 AVC, 5 HCC, and 
5 ICC. Samples were microdissected to maximize tumor cell 
percentage. In all samples an adequate library for sequencing was 
obtained. A mean coverage of 1800x was achieved, with 87.4% 
target bases covered more than lOOx and a mean read length of 78 
base pairs. 

The spectrum of mutated genes detected in the various tumor 
types was consistent with the current literature as reviewed in the 
COSMIC database [18]. KRAS mutations were detected in all 
PDAC and in 3 of 5 AVC and IPMN; TP53 mutations in 3 of 5 
PDAC, 3 of 5 GC and 2 of 5 IPMN, respectively. Other frequently 
mutated genes were GNAS in IPMN (4 of 5 samples), IDHl in 
ICC (2 of 5 samples) and PIK3CA in GC and HCC (Table 2). 

Twenty tumors (57%) showed multiple gene alterations with 
PDAC and IPMN displaying up to four concurring different 
alterations. In seven cases (20%), significant differences were 
observed in the frequencies of alterations affecting distinct genes, 
suggesting the presence of intra-tumor molecular heterogeneity. 
For example, case ICC4 had a 25% frequency Q61R mutation in 
the NBAS gene coexisting with a 9% oi BRAF Q461*, while 
ICC 1 showed a 18% frequency R282W TP53 mutation coexisting 
with a 5% oiSMAD4 C115Y. 

Orthogonal validation of intratumor molecular 
heterogeneity detected at TM-NGS by immunolabelling 
for p53 and Smad4 

To validate the relationship between mutation frequency and 
intratumor heterogeneity, we performed IHC analysis for p53 in 
ICC and for Smad4 in AVC samples. Case ICC5 (72% mutation 
frequency) showed a strong and diffuse immunostaining, whereas 
ICCl (18% mutation frequency) showed a heterogeneous pattern 



with sparse or clustered positive cells, roughly accounting for 15% 
of immunolabeUed cancer cells surrounded by regions of negative 
staining (Figure 3). As for Smad4 immunohistochemical analysis 
in AVC, the tumor sample AVC4 bearing a R361H mutation with 
24% frequency displayed a mixture of negative and positive 
regions, the latter accounting for about 15-20% of cancer cells, 
while non-mutated samples had a homogeneously positive 
immunohistochemical pattern (Figure 4). 

Discussion 

The results of our study may be summarized as follows: TM- 
NGS can be applied on DNA from routinely prepared paraffm 
tissues; the data produced are quantitative and thus permit the 
description of the molecular subclonal composition of a tumor. 

The introduction of targeted drugs is changing the profile of 
information needed to plan a therapeutical approach that entails 
multiple lines of intervention [19-22]. In this scenario, the 
histopathological diagnosis based on morphological classifications 
is no longer sufficient, and wiU need to be complemented by a 
comprehensive description of the specific molecular alterations 
and clonal heterogeneity of the tumor [23-25]. Proof of concept 
reports have already shown the potential application of NGS 
techniques using DNA from FFPE tissues [9,26-29] . However, its 
introduction in the clinical routine still needs validation of each 
step leading from the sample to results as well as the design of 
appropriate panels to specifically interrogate multiple tumor 
categories. 

The present study was therefore designed to evaluate the 
practicability of TM-NGS in detecting heterogeneity among 
diverse tumor types of the upper gastrointestinal system. In 
particular, three issues were addressed: i) to compare the 
performance of TM-NGS on FFPE-derived partially degraded 
DNA with that on high molecular weight DNA from fresh-frozen 
tissues; ii) to assess the mutation detection limit of TM-NGS on 
both fresh-frozen and FFPE derived DNA; iii) to assess TM-NGS 
ability in detecting inter-tumor and intra-tumor heterogeneity 
across upper gastrointestinal tract neoplasms. 

We used a commercially available multigene panel that 
simultaneously investigates the status of mutational hotspots of 
46 genes, including oncogenes with available (EGFR and BRAF) 
and upcoming [MET and PIK3CA) targeted therapies, or known 
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Figure 2. Sensitivity of TM-NGS for mutation assessment in fresKi-frozen and formalin-fixed paraffin embedded samples. A) DNA 

from AVCl cell line and a commercial germllne DNA were mixed to obtain samples with a decreasing relative AVC1 DNA content (50%, 25%, 20%, 
15%, 10%, 7.5%, 5%, 2.5%, 1% and 0%) to test TM-NGS sensitivity on DNA from frozen tissues. Three known homozygous variants harbored by the 
AVCl cell line {KRAS G12A, CTNNB1 S45F and the nonpathogenic polymorphism TP53 P72R) were used to assess the assay sensitivity. The commercial 
germllne DNA was heterozygous for the TP53 P72R nonpathogenic polymorphism. The variants were Identified In all samples containing AVCl DNA, 
down to a frequency of 1%. Obs: mutation frequency detected by Instrument, Exp: expected value calculated basing on dilution and mutation allelic 
frequency In the source AVCl DNA. B) A case of carcinosarcoma was used to test TM-NGS sensitivity In formalln-flxed paraffin embedded samples. 
DNA from two separate tumor components (hepatocarcinoma and sarcoma) was mixed to obtain samples with a decreasing relative 
hepatocarcinoma DNA content: 100%, 95%, 90%, 75%, 50%, 25%, 10% and 0%. These were subjected to the assay exploiting three known different 
mutations: ABLl Intronic g.164164 G>T (variant frequency hepatocarcinoma = 24% sarcoma = 0%), PIK3CA H1047R (variant frequency 
hepatocarcinoma = 28% sarcoma = 0%) and TP53 F109C (variant frequency hepatocarcinoma = 87%; sarcoma = 55%) The three mutations were 
Identified In all samples containing hepatocarcinoma DNA down to a frequency of 2%, corresponding to the frequency of ABLl gene mutation In the 
10% diluted sample. Obs: mutation frequency detected by Instrument, Exp: expected value calculated basing on hepatocarclnoma/sarcoma mixing 
ratio and mutation allelic frequency In each tumor component before dilution. 
dol:1 0.1 371 /journal.pone.01 04979.g002 



to decrease the efficacy of specific persoiiafized therapies [KRAS, 
NRAS, HRAS). 

Mutation detection by TM-NGS was as efficient with the 
partially degraded DNA from FFPE as it was with high molecular 



weight DNA from fresh-frozen samples, as shown by the similar 
coverage and allelic frequency of mutations obtained on matched 
samples of five SPT. The sensitivity of the assay was assessed by 
dilution curves, demonstrating that TM-NGS can detect mutated 
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Figure 3. The allelic frequency of mutation in TP53 gene corresponds to the proportion of p53 immunostained cells. Light blue 
sequence boxes indicate wild-type amplicons, light red indicates amplicons bearing a mutation, which is highlighted in red. Bars above amplicons 
show the relative abundance of wild-type (grey) and mutant (red) nucleotides. A) a case with wild type TP53 showing no p53 immunostaining; B) a 
case showing about 20% of immunolabelled cells for p53, consistent with TPS3 mutation frequency of 18%; C) a case with TP53 mutation frequency 
of 72%, showing a strong and diffuse p53 immunostaining. For each sample a representative H&E and p53 immunohistochemical image (original 
magnification x20) and the representation of the reads aligned to the reference genome are presented. 
doi:1 0.1 371 /journal.pone.01 04979.g003 



DNA accounting for 2 % of the cells in FFPE samples, reaching an 
even lower detection limit (1%) in fresh-frozen cells/tissues. 

While analyzing 35 samples from 7 different tumor types, seven 
cases with multiple mutations showed significant differences in the 
frequencies of alterations affecting distinct genes, while in ten cases 
the allelic frequency of mutations was not consistent with 
neoplastic cells percentage; this suggested the presence of intra- 
tumor molecular heterogeneity. Confirmation that TM-NGS 
quantifies the alleles affected, permitting the description of cancer 
subclonal composition was obtained by immunohistochemistry: 
this showed that p53 accumulation or Smad4 loss were seen in a 
proportion of cells comparable to that indicated by the allelic 
frequency of the mutation in the corresponding gene. 

The prevalence and type of mutations detected are comparable 
to those expected in the diverse tumor types considered herein, as 
reported by the curated COSMIC database [18]; the CTNNBl 
gene was always mutated in SPT and in 3 of 5 HCC [17], the 
R132 hotspot in IDHl gene was identified for ICC [30,31] and 
GNAS R201 for IPMN, KRAS was the most frequently mutated 
gene in pancreatic cancers while TP33 was frequently mutated in 
both pancreatic and gastric cancers [15,32,33]. Other frequently 
involved genes included PIK3CA and SMAD4. 



All the 35 tumor samples in our representative series of upper 
gastrointestinal system cancers were characterized by at least 
one single specific molecular alteration among the 46 genes 
analyzed, some of which also represent a potential therapeutic 
target. Two or more mutations were found in 20/35 (57%) 
cases. Moreover, several genes were altered in more than one 
tumor type, suggesting the possibility of a molecular subclassi- 
fication of tumors that crosses the borders of histology and puts 
the focus on molecular and potentially actionable alterations 
[34]. While these commonly altered genes could be detected by 
the commercial assay used in the present work, additional cross- 
border molecular alterations or mutations that remain confined 
to a specific tumor class are being reported [26,27]. For this 
reason, the design of specialized and optimized multigene panels 
will be the next mandatory step. Indeed, a European 
consortium of research centers has already developed a TM- 
NGS panel specifically tailored to target colon and lung cancer 
[14]. 

In conclusion, our study demonstrates the ability of TM-NGS to 
detect and quantitate multiple gene alterations, thus moving a 
further step towards a next-generation histopathologic diagnosis 
that integrates morphologic, immunophenotypic, and mutational 
analysis of multiple genes using routinely processed tissues. 
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Figure 4. SMAD4 mutational status corresponds to Smacl4 immunohistocKiemical loss of expression. Light blue sequence boxes indicate 
wild-type amplicons, light red indicates amplicons bearing a mutation, which is highlighted in red. Bars above amplicons show the relative 
abundance of wild-type (grey) and mutant (red) nucleotides. A) A case with wild type SMAD4 showing uniform Smad4 staining. B) A case with SMAD4 
mutation allelic frequency of 24% shows a heterogeneous pattern of immunostaining with alternating positive and negative areas. For each sample a 
representative H&E and Smad4 immunohistochemical image (original magnification x20) and the representation of the reads aligned to the reference 
genome are presented. 
doi:1 0.1 371/journal.pone.01 04979.g004 



Morphology and immunohistochemistry will provide diagnosis 
and drive the choice of areas to be microdissected for multiplex 
deep sequencing, while aiding the interpretation of sequencing 
data in light of intratumor heterogeneity. The role of the 
pathologist wiU be also critical to ensure the appropriate and 
ample sampling of the tumor to guarantee a complete and 
combined histopathologic molecular diagnosis. 

Finally, next generation targeted sequencing on paraffm tissue is 
much less expensive than the sum of many single conventional 
analyses, while having equal or even higher sensitivity [35-37]. 
This renders clinical application feasible and paves the way to a 
significant curtail of the economic burden of National Health 
Services. 
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